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Abstract 

A purely growing electrostatic drift instability driven by the electron tem- 
perature gradient at the inner edge of the plasma sheet can grow for large enough 
values of the temperature gradient. The parallel electric field associated with 
the instability is localized near the magnetic equator. The growth of the drift 
instability leads to enhanced whistler noise and increased electron pitch angle 
diffusion. If the current limit is exceeded in the ionosphere while the parallel 
electric field of the drift instability exists along the field line, rapid electron 
precipitation (the auroral breakup) can resvilt. 
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AN AURORAL BREAKUP MECHANISM 


This paper evaluates the possibility the auroral breakup is related to a 
purely growing drift instability at the inner edge of the plasma sheet. The 
auroral breakup here refers to the initial stage of an auroral substorm (T = 0), 
i.e, the sudden brightening at or near a pre-existing active auroral arc prior to 
the poleward expansion of auroral activity (see Akasofu (1968) for a discussion 
of the phases of a substorm). The inner edge of the plasma sheet refers to the 
near earth region of the electron plasma sheet where an exponential decrease in 
the electron energy density is observed (Vasyliunas, 1968; Schield and Frank, 
1970). 

Drift instabilities associated with the ring current particle population have 
been specifically discussed as causes of the auroral breakup by Swift (1967) and 
Liu (1970). The motivation for looking at a drift instability associated with the 
earthward termination of the plasma sheet as a cause of the auroral breakup is 
(1) its location on auroral field lines and its earthward movement with increas- 
ing magnetic activity (Frank, 1971), (2) the similarity between the energy spectra 
and energy densities of plasma sheet electrons and electrons precipitated in the 
auroral zone (Frank and Ackerson, 1971) and (3) the existence of a pressure 
gradient in the electrons at the inner edge of the plasma sheet. The purely 
growing drift instability is of interest because it becomes unstable at moderate 
pressure gradients and has an electric field component along the magnetic field. 


1 



Theory of the Drift Instability 

In discussing the time scales of the instability the following simplifying 
assumptions are made: (1) the geomagnetic field is axially symmetric, (2) the 
perturbation fields are electrostatic, (3) the ionospheres are perfectly conducting, 
(4) the temporal and spatial variations of the perturbation fields are large enough 
to preserve the first and second adiabatic invariants /u, and J. The limitations 
these assumptions place on the analysis will be discussed later, A method for 
studying drift instabilities in a general magnetic field configuration has been 
developed by Rutherford and Frieman (1968a) but with the assumption made 
above it is convenient to use a variational principle technique as used by Rosen- 
bluth (1968) or Liu (1969a), The analysis of the instability will closely parallel 
that of Liu (1969b). 

Euler coordinates a and /? will be used to describe the magnetic field 

B(r, (9, 0) = VaxV/3 (1) 

with an axially symmetric field it is judicious to choose ft to be the azimuthal 
angle, 0 , so that a. - a(r,9). The third coordinate of the system, denoted by the 
parameter s, is the distance along a field line from a fixed reference surface 
through which all field lines pass. 

In discussing the time scales of the drift instability the pertinent frequencies 
are the gyrofrequency, the bounce frequency, the azimuthal drift frequency and 
the azimuthal diamagnetic drift frequency. The gyrofrequency is 

Q. 

^ m. c 
i 
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where B is the magnetic field strength, e is the particle charge, c the speed of 
light, and m 3 the particle mass. The bounce frequency of a particle between 
mirror points in the magnetic field is 

_ BH(g, /3, M. J) _ r 

where 

H = K + e$ + — ^ (2) 

C dt 


and K is the particle kinetic energy and <D the electric potential. H does not de- 
pend on the distance along the field line and for bounced averaged guiding center 
motion cH/e plays the role of the Hamiltonian (Northrop, 1963) 


<a> = - 


c 

e 3/S’ 





The brackets < > denote a time average over the bounce period of the particle 

/V r ds r vox 

<fx> = V. I — ^ ^ 

Ms v,,(s) 


The particle velocity along the field line is denoted by vj| (s). The limits of the 
integration over s are the mirror points s_^ where v^^ (s^ ) = 0 and the magnetic 
field has the value B(s^) = B^. It is convenient to define 1/ v^| (s) to be zero for 
those values of s along the field line that the particle does not traverse in a 
bounce period, that is, where B(s) > B^. Then the integration in equation 3 can 
be considered to be over the total length of the field line. The azimuthal drift 
frequency is 
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</ 5 ) = 


"h = - S 

e da 


and the azimuthal diamagnetic drift frequency is 


o> (4) 

e Ba 

where p is the plasma pressure, T the plasma temperature and k denotes 
Boltzmann’s constant. 

In the magnetosphere near the inner edge of the plasma sheet the ion gyro- 
frequency in the equatorial plane is on the order of 1 Hz. The bounce frequency 
for a 5 Kev proton mirroring near the ionospheres is on the order of 10“^ Hz. 

For a 5 Kev particle at 6 the azimuthal drift frequency is on the order of 
10 ® Hz. The diamagnetic drift frequency depends upon the pressure gradient 
and is here considered to be on the order of the m^nitude of the azimuthal drift 
frequency. 

The instability associated with the inner edge of the plasma sheet is analyzed 
under the assumptions it is electrostatic and preserves /j. and J. The wave length 
of the perturbation field perpendicular to the geomagnetic field, Lj^ /2n , and par- 
allel to the geomagnetic field, L|| /27 t , as well as the scale height of the pressure, 
Lp , are taken to be of the same order of magnitude 

L -L,| :.L 

The fields and frequencies associated with the perturbation are ordered with 
respect to the parameter e = p/L, p being the ion gyroradius. The frequencies 
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are taken to be ordered as 


a> ~ e^Q. 

bi 1 

where w is the wave frequency, the ion bounce frequency and a the ion 
gyrofrequency. The fields are ordered as 

E - E, - e — (5) 

where E ^ and Ej^ are, respectively, the components of the perturbation fields 
parallel and perpendicular to the geomagnetic fields. The requirement on the 
perpendicular electric field states that the electric force is of higher order than 
the Lorentz force. The requirement on the parallel electric field means AV|| , 
the change in v,| due to E|| in one gyroperiod, is of higher order than 
V|| (A V|| = (e/m) E|j /fi = cE|| /B = fc V|| ). Note the parallel and perpendicular 
fields are assumed to be of the same order so the geometry appropriate to these 
waves is not the long thin geometry of the finite Lamor radius ordering in which 

E|j e2 2 Lb; Lj^~eL||i 

The existence of an electrostatic instability can be determined by investi- 
gating solutions to the coupled Poisson and Vlasov equations. In the ordering 
presented above, Poisson’s equation is equivalent to the charge neutrality con- 
dition to second order in e. Poisson's equation is 

V- (E|| + Ep 47re(n. - n^) 

where n^ is the ion number density and n^ is the electron number density. With 


5 



the fields ordered as in equation 5 


Bn 


V-E 


and 


Bn 


477-e(n. - n^) 


c 


1 




L|i 

+ — 
L 




_ 1! 

nj 



'p " 

2 (n 

i “ '"e 

:) nj - 

"l “ 



n. 

1 

n . 
1 



where — [kT/4?rnje^] is the Debye length. Then Poisson's equation requires 


n. - n /t \2 

For temperatures, densities and field strengths near the inner edge of the plasma 
sheet 

and n. - n < < e^n. 

Thus to order e ^ Poisson's equation can be replaced by the charge neutrality 
condition 


^(r>i - "e) = 0 

When the first and second adiabatic invariants are conserved Vlasov’s 
equation can be written (Northrop, 1963) 


^ + ii f— ^ M 

dt e 9/3 9a 9a 9/? 


( 6 ) 


where F - F(u. ,J,a,yS,t) is the phase space distribution density. Denoting the 
equilibrium quantities with a subscript "O” and perturbation quantities with a 
subscript ”1” F can be written 


F(p, J, a, A t) = F^(p, J, H) + Fi(/a, J, a, /S, t) 
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The quantity (m ,J,H) can easily be seen to be an equilibrium solution (3F/^t= 0) 
to equation 6. With axial symmetry it is convenient to write the perturbation 
quantities 0^, • Fj in component form 




s) exp i(Z/3 - o)t) 


I 


Hj=^^H^(a, s) exp - £<jt) 

I 

Fj = ^ _' fj) (g, s) exp i(-tj3 - <nt) 


( 7 ) 


where a time dependence exp -iwt has been assumed for the perturbation quan- 
tities. The condition J is conserved is interpreted in the sense <j) = 0. Using 
J conservation and noting that with d/3/dt ~ 0 


J = <b dsv.| (s) = d) ds 


- (H - - e$ ) 

.m 


1/2 


the components of Hj can be related to the components of 


h = ^ j ^ s) 

The spatial number density of particles n(r,t) at the point r = , s') is related 

to the phase space distribution density as (Northrop and Teller, 1960) 

n(7. „ 

= 2B [d,ifdJv;'(M. J. s') F(/x, J, a', /3', t).-,^ (9) 
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The integration over , J in equation 9 is restricted to those values of /u., J for 
which the bounce motion of a particle carries it through the point (a' , /3' , s') = r. 
However, with the extended definition of v'^ (s) discussed above the integration 
will be taken over all/u. , J since v”* (;u, J, s ) = 0 for particles with ^ , J values 
such that they do not pass through s'. 

Using equations 6 and 7 the components of the perturbed distribution density, 
Fj, can be expressed in terms of the components of the electric potential, 
and the equilibrium distribution density (yLi,J,H) 







M.j 



or writing a> =co' + iy with o)' , y real 





Using equation 9 and 


( 10 ) 


the components of the perturbed number density can be written 



Since - e and H is independent of s, is a function of H, a. only, 

and (after Rutherford and Frieman, 1968b) 
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or performing the indicated differentiation and using 


^ =0 and - v-3 = tn — (v-i) 


3H 


" 3H 


fi , a 


- 


I <$£>* •'V li * I — 3 

^ mv,r ^ ^ mv^ 


= 0 


or 




mv3 

Using the above equation and = 9 J/9H equation 11 is 

= 2B|d^JdJ 


/X , a 


F^- e 


0 3j 3j 


or integrating by parts 


no r 2B Id^i 


r 


fi, a 


The charge neutrality condition 

j 

(where S is used to denote a sum of the electron and ion terms) along with 
j 

equation 10 gives 
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The ionospheres are assumed to be perfectly conducting so that the boundary 
conditions associated with equation 12 are and d^/ds both equal zero at the 
ionospheres. If <I>^ is constant at the ionosphere, dO^/ds must be zero since 
the field lines are not normal to the ionospheres. 

The ordering used to analyze the Vlasov equation has resulted in an eigen- 
value equation for localized in a (eq. 12). A variational principle can be 
formed from equation 12 and used to obtain information concerning the eigen- 
values and eigenfunctions by multiplying equation 12 by the adjoint solution, , 
and integrating over the length of the field line (see Morse and Feshbach, 1953). 
The adjoint solution to equation 12 can be foimd by investigating the nature of the 
kernel of the integral operator C , where eq. 12 is written 

£<I)^ =: and -J dsL(s'/s) $^(s) 

The limits of the integration over s are now fixed at the ionospheric values of s, 
and the kernel L is 


L(s'/s) 



dj g(M. J) 


J) 

v„Cs') 


J) 

V||(S) 


wdiere g(A<-, J) does not depend on s and 
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The operator C is self-adjoint since L(sl s') - L(s '1 s) but it is not Hermetian 
since L*(s| s') ?^L(s'l s) where denotes the operation of complex conjugation. 
Since the boundary conditions are adjoint and L(s] s') = L(s'l s) the adjoint equa- 

tion is the same as equation 12 and ij)^ = 

The variational principle which yields equation 12 is 8D/8$^ = 0 where 



The eigenfunctions, <I>|, extremize D, i.e. 5 D/S^>^ consistent 

with equation 12 the frequency must be such that <i>|) = 0. It will be shown 
below that under some conditions there exists an u> with a positive imaginary 
part (7 > 0) and a<J)| such that (D(w ,<I>|) = 0 and D is stationary at <i>l . Thus 
is an eigenfunction of equation 12 with a corresponding w such that it is purely 
growing in time. 
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Conditions for the Existence of the Purely Growing Drift Instability 

For the purpose of illustrating the physical processes contributii^ to the 
growth of the instability and simplifying calculations the magnetospheric plasma 
is assumed to have an isotropic pitch angle distribution and a Maxwellian dis- 
tribution in energy. A more realistic evaluation of the stability criteria using 
experimentally determined differential electron energy spectra will be made 
later. To express the equatorial number density it is convenient to make a change 
of variables from J to A. , K where X. = /x/K is the inverse of the mirror point 
magnetic field strength. The equilibrium equatorial number density n^ (a) can 
be expressed as an integral over A , K 


n.(a) - 2B^ 



dj 


- 

F^(m, J. a) = 



X ^(1 ^ 2Bg]/5 


dAdKK^/2 
(1 - 


K, a) 


where with the above assumptions and noting that the charge neutrality condition 
implies n^(a) = n.(a) - 


F-(K. K, a) = "oWexp-[K/kT,(a)] 


B;' £ \ 

E -s J. 


(277m.k3Tj>(a))l/2 

k is Boltzmann's constant, T. (a) is the temperature of the jth species, B^/Bj 
is the ratio of the equatorial to ionospheric magnetic field strength which 
taken to be small enoi^h to write 


IS 
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F^J(K, a) = exp - [K/kT (a)] 

(277mj k^T?(a))^^^ 

Also, using the variables ,K), the bounce frequency and the azimuthal drift 
frequency are factorable in terms of the particle kinetic energy, K. If the 
equilibrium electric field is zero 




2K 

m . 

J 


|l/2 




K 


kT. (a) 


is the azimuthal drift frequency of a particle of the jth species with energy 


kT. and 


1/ - f ^ 

is independent of the particle species. 

Equation 12 will be examined for eigenfunctions corresponding to purely 
growing modes (that is oj' =0). It will be shown that a class of functions exist 
that extremize D where the extreme values of D can satisfy D = 0 under some 
conditions with y > 0(o) = iy). The azimuthal drift frequency co^ will be taken 
to be independent of \ as is the case for a dipole field. Equation 12 for the 
purely growing modes and a Maxwellian energy distribution ) is 


D - ^gj(a, y2) + ijgjCa. 7^) - gaCa)^ + 



where 


giCa. 7^) 







"d 


2 


CO'- + 


c 3K 


3FA 


+- 


/i. a 


M.Kj 


^2 
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r 1/2 


Q 


o-^l + 

^ da 

" 



c BK 




1- 

/X, 01 





-e 2 


g,C 





K>/ 2 dKe. ~1 

^ 3K 


M> ^ 


The functions g. are all real. 


Now let = R + i I where R and I are real functions of s, so that 

<♦[>' =<RV -<t>^ t2i<R><}> 

<®|> = <E‘^-<I2> + 2i <RI> 

and with D = + i Dj {D. , D. real) 


D = I d\ V- 1 (^) « R>' -<'r>2 ) (g, - g3> - 2 I g2<R><!> + g3 «R 2> - <I 2 >)] 


and 


D. = Jd\z.-i(\) |'«r> 2 . <f>2 ) r ^ 2<R><I>(g, - 63) + 2g3<RI'>] 

It can now be shown that there exists a class of functions that simultaneously 
extremizes and D^. The functions which accomplish this are those with R = al 
where a is a real constant, since then 



a^) (gi - gg) - 2a ^ gj + g3<R2>\ 


and 



(1 - Jga + 2 a< 1 ^^ 


(gj - gg) + 2a 
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Then the functions which extremize both and D. simultaneously are those 
which extremize the quantity 

b = 

|d\v-i<Ry 

Let the extreme value of b equal M then the conditions - D. = 0 for D station- 
ary in are 


- 0 - (1 - _ 2^g^a 


(14a) 


D. ^ 0 = (1 - + 2agj + 2agj(M - 1) (14b) 


which are two equations in the real variables a and y. 

Solutions to equation 14 can be easily found for marginal stability (7 = 0) 
provided that the plasma pressure gradient satisfies a certain condition. Choos- 
ii^ a = 0 (i.e. real) gives D. = 0 and D = 0 if 


gg(a) CM - 1) + gj(a, 72 = 0) = 0 

or with the definitions of g^ and g^ above the condition for marginal stability is 





d InrijjCa) 
da 



M 

2 


(15) 


J 

It can be shown that solutions exist for y ^ 0. Suppose 7 is small and ap- 


proximate g^ by 


gi(r^) 


g/O) 


1 +2L 

1^2 


-= gjCO) - 2 2_ g (0) 
2 ^^2 
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and note gj is independent of y. Equation 14a and 14b are 


(1 a^) g 3 (M- 1) + (1 - a^) gj(0) - 2 — (i _ a^) g (0) + 2ZgJy) a = 0 

^2 ^ I 


(1 - a^)^g^(y) + 2a[gj(0) + ggCM - 1)] - 4a21 g^(Q) = 0 


the first equation then gives 


(1 - (gaCM - 1) + g,(0)) - 2 2l!gj(0) (1 - a2) /2 


an expression that can be satisfied by either sign of y simply by changing the 
sign of 'I . Using this result in the second equation gives 


[gjCM - 1) + gj(0)] 

i/O) 


(1 . ay - 2(1 - a^)^2i: . 8aV = 0 


SO that if gj{0) > 0 solutions exist if g (M-1) + g (0) > 0 or 

3 1 


jy-c f -dlnya) dlnTI I „ 

kT. [da da ^ ^ ^ [ 9 


It has been shown that the class of functions f — (1 +ia) a real function, 

a real) can simultaneously extremize the real and imaginary parts of D and for 
certain values of the plasma pressure gradient the condition D = 0 can be met 
and is an extremum for the function f’ = (1 +ia^) which extremizes 

< 4 .|> 
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where satisfies the boundary conditions noted above d4>^/ds = 0. Accord- 
ing to the Schwartz inequality (with1>|^ real, = l^| I ) the mininaum value 

of b occurs for > a value achieved if is a constant along the field 

line. However a non-zero constant potential is ruled out by the boundary condi- 
tions. So the minimum value of b will be somewhat greater than unity. 

The fundamental or fastest growing eigenfunction for each value of 't is a 
symmetric function of s (i.e. centered about the magnetic equator s = s^) and 
zero at the ionosphereic values of s. The dependence on s of the magnitude of 
the fundamental eigenfunctions is determined by the requirement b is minimum. 
Consider a potential fmiction O' with scale size L'. The average value of O' can 
be approximated for a given k. by 


<0'> ~ f(\) 0^ ^ 


L' 




where s^ is a value of s near the turning point so that f(X) can be interpreted 
as the fraction of the bounce period spent within a scale size, L’, of the turning 
point. The factor l/V|| in the integrand of the averaging operation causes the 
potential in the vicinity of the turning point to make the largest contribution to the 
average value. From the Schwartz inequality<C|^| ^ ^ S )>| ^ it follows 
f(^) 5 If for all>^ , and then 



d\y-if(^) O^/f 0^] 


(16) 


is smallest when f(\) is unity. f(\) is only unity for particles mirrorii^ at the 
magnetic equator, b will be minimum for potentials with magnitudes maximum 
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at the equator and falling rapidly to zero since then f(^) will be near unity when 
the potential is non-zero in the integration in equation 16. So the fastest growing 
potential eigenfunctions will be concentrated near the magnetic equator falling to 
zero in a distance on the order of the equilibrium pressure scale height. 

Using a dipole field to approximate the geomagnetic field = (-3/a) (c/e.)kT. 
and assuming the temperatures of the species are equal the condition for the 
existence of a purely growing instability becomes 


dlnng(a) dlnT(a)1^3., 3 

dlna dlnaj 2 2 


(17) 


Gradients in the temperature and number density determine the stability of the 
equilibrium. At the inner edge of the plasma sheet the temperature gradient will 
be found to be the most important factor in determining the stability of the plasma. 

The electric fields associated with a purely growing instability will mostly 
affect the orbits of low energy particles (o)^< wJ). The instability derives its 
energy from particles which are moved about in the geomagnetic field by the 
perturbation electric fields of the instability. If the net energy of the particles 
interchanged by the wave fields is less than before the interchange, energy is 
available for the instability to grow. In the familiar hydromagnetic interchange, 
instability onsets if a sufficiently large gradient in the number density exists in 
the same direction as the magnetic field gradient. For the purely growing in- 
stability equation 15 indicates the temperature gradient must be opposite the 
field gradient for instability to onset. This behavior of the temperature gradient 
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arises from the term 


rdKKi /2 r 3] dlnT 

Jw^K/kTLkT 2J da 

in the evaluation of D at marginal stability. It might be thought that the sign of 
the gradient in the temperature indicated in equation 15 as leading to instability 
should actually be stable due to the decrease in the average energy of the par- 
ticles with decreasing temperature. Actually the temperature gradient leads to 
instability in a manner quite similar to the familiar hydromagnetic interchange 
in that the number of interchanged particles (in this case the low energy par- 
ticles) increases for increasing B. The low energy particles increase due to the 
change in the normalization factor l/T^/^ with a (this chaise is expressed by 
the term 



in the evaluation of D). The change in the normalization factor more than com- 


pensates for the decrease in the average energy of the particles with decreasing 


a (the decrease in average particle energy is represented by the term 


dlnTfdK K K 

da J jji/2 kT kT 

in the evaluation of D). As a result the energy change favors the growth of the 
instability for 
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(in a dipole field) and the low energy particles in the plasma supply the free 
energy for the growth of the instability. It might be expected, then, that non- 
Maxwellian distributions in the energy of the plasma \^ich have a similar be- 
havior (i.e an increase in low energy particles with decreasing a ) may also be 
subject to the purely growing instability. 

The requirements for the onset of the instability can be expressed in terms 
of the directional differential intensity, I (a quantity commonly used to present 
satellite data). The directional differential intensity is related to the phase 
space distribution density, by « I/K. Using this expression in equation 13 
leads to a requirement for marginal stability in a dipole field of 




fdKK-3/2 I 

J l3a '^kJ 


= M ~ 1 



(18) 


The expressions involving I in equation 18 can be approximated by the directional 
intensity measured for several energy bandpasses as a fimction of radial distance. 
Such measurements for electrons at the inner edge of the plasma sheet are avail- 
able from satellite experiments (see e.g. Schield and Frank, 1970; Frank, 1971). 

The growth rate of the instability with F Maxwellian can be found by solving 
equation 13 for 7 . Unfortvinately the integration over the energy is not trivial. 

The growth rate can be estimated by approximating the integrals in equation 13 
by, for example, 
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the last estimate foUowing since the diamagnetic drift frequency is on the order 
of the azimuthal drift frequency. Note the growth rate does not maximize be- 
cause particle losses arising from violation of J have not been accounted for. 
The assumption of m > J conservation requires that the distance a particle drifts 
a^imuthaHy in a bounce period is small compared to the perpendicular wave 
length of the pertiurbation fields 


rw < <K7i « 

Ob i ^ 

so that . BvLt so that even for short azimuthal wave lengths the 

0 b 0 

fastest growth rate to be expected is on the order of 10 sec . The growth 
rate of the instability in relation to the auroral breakup will be discussed later. 

In evaluating the proposed instability as a cause of the auroral breakup it is 
important to consider the factors which tend to stabilize the plasma. Objections 
have been raised to ascribing auroral breakup to electrostatic interchange 
instabilities (e.g. see Hasegawa, 1971) since they are easily stabilized by a small 
amount of homogeneous cold plasma. If it is assumed that such a plasma is 
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present at the inner edge of the plasma sheet with number density n^ and tem- 
perature then the requirement for instability with f Maxwellian becomes 


din. 2 

If a cold plasma is to stabilize the interchange it should be noted it is neces- 
sary that the cold plasma be homogeneous as illustrated by the following example. 
Suppose the inner edge of the plasma sheet is formed when plasma sheet electrons 
are precipitated into the atmosphere and replaced by cold ionospheric electrons 
(somewhat heated during the replacement process) in such a way that the inner 
edge is a two component electron plasma with Maxwellian velocity distributions 
and with number densities and temperatures n , T and n,T where T <<T. The 
electron number density at the inner edge is constant, n ^ + n = constant. The 
apparent decrease in temperature at the inner edge is now due to an increase in 
n^. In this case 

da da 

and the stability criteria is 

\T / dlna 2 
At the inner edge of the plasma sheet 



and the two component plasma could be unstable. The addition of a cold electron 
component to the plasma from the ionosphere might take place in such a way as 
to lead to a situation in which the added cold plasma eventually becomes the 
source of free energy for a purely growing instability. The existence of a cold 
plasma component in the vicinity of the inner edge of the plasma sheet has not 
been established by direct satellite measurement due to the experimental diffi- 
culties involved. In the following a homogeneous cold plasma component is 
assumed absent at the inner edge of the electron plasma sheet. 

An anisoptropic pitch angle distribution also alters the stability criteria. 

For a fixed temperature gradient a flat pitch angle distribution (i.e. a distribu- 
tion having more particles with velocities perpendicular to the magnetic field than 
an isotropic distribution) is more stable than a field aligned pitch angle distribu- 
tion. The alteration of stability criteria is due to the energy durir^ an interchange 
of particles at a fixed temperature gradient being larger for averse larger J. 

The effect is small for pitch angle distributions not too far from isotropy and is 
not considered here. 

The effect of an equilibrium electric field with potential P(a) depending only 
on CL can be established from equation 12. Such a field is representative of the 
corotational field of the earth or convection-like fields in the dawn and dusk 
sectors. The effect of such a field is to make the equilibrium overstable. That 
is, the fastest growing mode corresponds to an eigenvalue with the real part of 
the frequency 
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, _ dP _ c dP 
da RB dR’ 

where K is the equatorial radial distance corresponding to a. A non-local analy- 
sis in a would introduce terms dependent upon d^ P/d^ a and are not considered 
here. The physical process leading to the growth of the instability can be learned 
by investigating the zero equilibrium electric field case. 

Stability of the Inner Edge of the Plasma Sheet 

The electron plasma sheet is a population of electrons in the magnetosphere 
of average energy 1-5 keV forming a layer some 4-6 R in thickness centered 

e 

about the neutral sheet (Vasyliunas, 1968; Schield and Frank, 1970). It extends 
near the earth into the dawn and dusk sectors where its thickness increases. 

The inner edge of the plasma sheet is in the portion of the magnetosphere where 
the magnetic field is nearly dipolar (Frank, 1971). The termination of the plasma 
sheet is identified by a sharp drop in the energy density of the electrons. The 
drop in the energy density is due to a decrease in the number of higher energy 
particles not a decrease in the number density. The plasma sheet terminates in 
the vicinity of 6-10 R^ from the earth and the location of the inner edge depends 
on Dst. The inner e(^e moves earthward as the level of magnetic activity in- 
creases. The plasma sheet electrons are characterized by their energy spec- 
trum which shows a peak intensity for energies between 1-5 keV. The energy 
spectrum characteristic of the inner edge of the plasma sheet is similar but has 
a break in the slope of the higher energy portion of the spectrum. The change in 
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the shape of the plasma sheet energy spectrum at the inner edge is shown in 
figure 1 (Figs. 4 and 7 of Schield and Frank, 1970). Nearer the earth the niunber 
density of particles, both protons and electrons, in the energy range 1-50 keV 
drops off at what is usually identified as the inner edge of the ring current. 

Usually the decrease in electron energy density occurs in the region where the 
number density is changing slowly. Thus the drop in energy density can be thought 
of as a decrease in the electron temperature. The important contribution to in- 
stability can be considered to come from the temperature gradient in the elec- 
trons. The scale length associated with the drop in electron temperature is 
.6-.4 R (Frank, 1971; Vasyliunas, 1968) and it occurs between the trapping 
boundary and the plasmapause. Some examples of the location of the decrease in 
plasma sheet electron energy density relative to the plasmapause, trapping 
boimdary and proton ring current are shown in figure 2 (Figs. 1 and 5 of Frank, 
1971; note also his Figs. 2-4,6). 

The directional differential intensity of electrons for various energy band- 
passes as functions of radial distance near the equator is shown in figure 3 (fig- 
ure 2 of Schield and Frank, 1970). Note the increase of low energy particles with 
decreasing L value (decreasing a) as well as a notable lack of particles with 
energies below 500 eV compared with a Maxwellian distribution (see figure 4). 

The use of equation 17 with Maxwellian energy distributions to represent the data 
of figure 3 as shown in figure 4 yields a value of approximately 10 for the left 
hand side of equation 17. This method would predict the inner edge was unstable 
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but the Maxwellian approximation is not a good choice due to the absence of low 
energy particles in the measured spectra. The departure of the distribution from 
Maxwellian makes it desirable to use equation 18 in evaluating the stability of the 
inner edge. A numerical calculation of the left hand side of equation 18 using the 
detailed data shown in figure 3 yields a maximum value of approximately 1 at 
6.7 The data in figure 3 corresponds to figure 2a. and for this particular 
case the electron temperatvire gradient is near the inner edge of the ring current 
and there is probably a stabilizing contribution from the protons. The tempera- 
ture drop in electrons does not usually occvir in the region where the proton 
number density decreases so figure 2a. is not a "typical" example. The impor- 
tant thing here is to note that a realistic estimate of the stability criteria shows 
the inner edge is near marginal stability. It is not unreasonable to expect that 
the increase in low energy particles generally observed at the inner edge of the 
plasma sheet can be large enough that the purely growing drift instability de- 
scribed above can be expected to grow and eventually result in the relaxation of 
the temperature gradient. 

Precipitation of Particles 

The growth of the instability alters the distribution of particles in phase 
space. A distribution of charge (less than a-t any point along the field line) 
sets up a parallel electric field. As discussed above the potential associated with 
the fastest growing mode will be non-zero only in the region of the magnetic 
equator. The change in the equilibrium phase space distribution density resulting 
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from the growth of the instability depends upon the values of \ , K and s. The 
equation giving the change in the spatial number density at any point along the 
field line is 


n^(s) = 2 B(s) 


rr dXdKK^/2 J dFg 



2jj(l-\B(s))i/2 I da 
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Note that replacing the term contained in braces in the integrand of equation 19 
with the equilibrium distribution yields n^(a) (a constant along the field line). 
Thus to compare the relative magnitudes of the changes in the equilibrium num- 
ber density along the field line it is sufficient to compare the magnitude of the 
terms in braces in equation 19. The largest changes occur in the vicinity of the 
magnetic equator where and are large. Changes in the velocity space 
distribution for values of s near the ionosphere are a factor 
{"^e) less than the changes in the equilibrium distribution near the equator. The 
reason for the small change in the near ionospheric velocity space distribution is 
that the parallel electric field of the instability grows slowly enough that the first 
and second adiabatic invariants are conserved. As the instability grows the 
energy of the particle changes so that its J is conserved. Denoting equilibrium 
- 0) quantities with the subscript "0," J at ^>^ 7 ^ 0 is 
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Since J does not change it is required that 


The change in the pitch angle, 6 = tan"* (Vj^ /V|| ), of a particle mirroring near 
the ionosphere with energy kT will then be 



Thus even if the change in the potential near the equator is on the order of kT 
few particles are dumped into the loss cone, since the pitch angle of particles 
mirroring at the ionosphere changes from7r/2 to tt/2 + [ (ci)^)/k] . Actually 

the pitch angle distribution of electrons and protons in the region along the field 
line where the parallel electric field is zero must be the same. If it were not 
the pitch angle distribution would then cause a non-zero electric field along the 
field line in that region. So any particles dumped into the loss cone consist of 
equal numbers of protons and electrons. The number of particles dumped, 
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however, is small. The instability can be considered to grow, setting up a paral- 
lel electric field near the magnetic equator, without precipitating particles. 

An upper limit to the magnitude of the- parallel electric field associated with 
the instability can be found from the ordering requirement E|| ~ e (vB/c). This 
requirement assures that the change in the parallel kinetic energy due to the 
electric field applied during one gyroperiod is of order ekT. An electric field 
consistent with this ordering could cause a change in the parallel kinetic energy 
of a particle 


AW„ = eA<E>= eE,,L= efl-kT 

but L/v - and eQ so that AW|| ~kT. The electric field of the instability 
can be large enough to change the parallel kinetic energy of particles a substantial 
amount as they move along the field line in the region of the magnetic equator. 
However, as discussed above, a potential change even of order kT does not cause 
rapid precipitation of particles. 

Although the instability does not immediately precipitate particles the inter- 
action of electrons with whistler-mode noise can be expected to diffuse electrons 
into the loss cone (Kennel and Petschek, 1966). The growth of the purely growing 
drift instability modifies the equilibrium in such a way that whistler-mode noise 
should be generated. The velocity space distribution will be unstable to waves 
of frequency and wave number k' if the anisotropy in the electron distribution, 
A" Petschek, 1966, where 
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F~ is the electron velocity space distribution and 6 is the pitch angle of the 
particle {&= tan"^ (v^ /v^ )). The change in the equilibrium equatorial number 
density is given approximately by {equation 19 with 72/^2 ^ 2 ^ 


"e-e = 


dX 


(1 - 




?<<!>> 


d In 


a kT V 3 d In a 7 


( 20 ) 


L 2K d In , 

The electron velocity space distribution is related to the term in brackets in 
equation 20. The distribution in velocity space after the instability has started 
to grow will depend on V (since <<i)> depends onX ) and thus 0 since = sirf 0 . 
Since the effect of the growth of the interchange is to relax the temperature 
gradient and decrease the number of low energy particles when 
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choose > 0 for low energies (K < kT). For this to be the case it is required 
that <^)> > 0. The magnitude of for the fastest growing mode decreases 
monotonically from large X (particles mirroring near the equator) to small X 
(particles mirroring near the ionosphere). Since <( 4 »> > 0 , <' 1 >> is monotonically 
decreasing as X decreases (and thus 0 decreases) so that 
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Thus for parallel energies greater than 3/2 kT, BF/^0 is positive for all K and 
A“ > 0. For high enough resonant energies. Eg, the condition for instability be- 
comes approximately, A“> 0, and so some whistler noise should be generated at 
frequencies 

where E_ is large enough that A~ > w' /(cj' - fi ) where the frequency satisfies the 
above expression and E,, = B^/87rn is the m^netic energy per particle. A", for 
a given E , depends on the magnitude of the potential associated with the par- 
allel electric field of the drift instability. The larger the potential the wider the 
range of parallel energies which satisfy the criteria for instability and the larger 
the growth rate. Thus whistler noise should be generated over a wider band 
width and interact with more particles as O increases. For the high energy 
portion of the measured plasma sheet spectrum the Maxwellian is a good approxi- 
mation and the above arguments can be expected to hold. 

If whistler noise is present in the vicinity of the equator it can be expected 
to diffuse some electrons into the loss cone. If a parallel electric field is present 
the loss of these electrons will change the potential in the vicinity of the iono- 
sphere. The reaction of ionospheric electrons to the changing potential near the 
ionosphere will be important to continued maintenance of an electric field along 
the field line. The loss of some electrons to the atmosphere will cause a tem- 
porary excess of positive charge along the field line near the ionosphere. This 
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charge will be neutralized by the highly mobile ionospheric electrons thus main- 
taining zero electric field in the vicinity of the ionosphere „ It will be assumed 
here that as long as ionospheric electrons maintain their mobility a small amount 
of plasma sheet electrons can be lost by pitch angle diffusion without altering the 
condition of zero electric field near the ionosphere. As long as this particle loss 
is small the purely growing drift instability should continue to grow with the same 
potential shape, i.e. non-zero only near the magnetic equator, described above. 

If the diffusion of particles into the loss cone is large and the instability has 
grown to a point where the parallel electric field is large, rapid precipitation of 
electrons may occur. If the number of electrons diffused into the loss cone by 
whistler noise is large enough that the current of neutralizing ionospheric elec- 
trons exceeds a certain limit, a two stream instability will onset in the ionosphere 
(Swift, 1965; Kindel and Kennel, (1971). The dominant modes growing as a result 
of the current limit being exceeded depends upon ionospheric composition and the 
magnitude of the current as has been shovm by Kindel and Kennel (1971). When 
the current is exceeded and instability onsets the resultii^ wave turbulence 
along the field line greatly impedes the ionospheric electron mobility. Under 
these conditions the ionospheric electrons cannot completely neutralize the 
excess positive charge near the ionosphere which results when plasma sheet 
electrons are precipitated into the atmosphere. As a result electrons c mning 
down the field line will "see" more positive charge near the ionosphere than 
would exist if no electrons had been precipitated. As a result their parallel 
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velocity increases and some are dumped into the loss cone. Since electrons 
move much faster than ions along the field line the ions can be considered fixed 
during this process. Dumping of electrons into the loss cone sustains the cur- 
rent and wave turbulence as well as increasing the excess of positive charge in 
the vicinity of the ionosphere which in turn leads to more electron precipitation 
and so forth, A large number of electrons will be rapidly precipitated until the 
ions redistribute themselves along the field line. This redistribution can be 
expected to take at least a few ion bounce periods (on the order of 5 minutes). 
The rapid precipitation which occurs if the critical current limit is exceeded 
when the parallel electric field of the purely growing drift instability exists 
along the field line is proposed here to be the auroral breakup. 

Discussion 

The auroral substorm is a part of the magnetospheric substorm which 
occurs in the polar ionosphere during magnetic storms and lasts for 1-3 hours 
(Akasofu, 1968). The magnetic perturbation associated with the auroral sub- 
storm and the latitudinal and longitudinal extent of the auroral substorm can 
vary greatly. The magnetospheric substorm has been proposed to consist of 
three phases; a growth phase (McPherron, 1970); an expansion phase; and a 
recovery phase (Akasofu, 1964). The growth phase has been cited as a period 
of enhanced magnetospheric convection during which (among other phenomena) 
the inner edge of the plasma sheet moves earthward, the plasma sheet thins, 
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and a current system similar to the 3-dimensional electrojet current system 
(Akasofu and Meng, 1969; Bonnevier, et, al,, 1970; Kisabeth 1972) is established 
(McPherron, 1972). The evidence presented for the existence of a growth phase 
has been disputed by Akasofu and Snyder (1972). The expansion phase marks the 
beginning of the auroral substorm as described by Akasofu (1964). Three main 
events are assigned to the expsinsion phase; (1) the sudden brightening near mid- 
night at or near a pre-existing auroral arc (duration 0-5 minutes); (2) the auroral 
arcs become distorted and move poleward (5-10 mins.}; (3) the westward travel- 
ing surge and auroral bulge is formed with the northern-most arc the most in- 
tense (10-30 mins.). The recovery phase begins when the arcs reach their most 
poleward position and start to recede equatorward and is marked by pulsating 
aurora and eastward drifting patches (1-2 hours). The rapid precipitation of 
electrons discussed in this paper deals with event 1 of the expansion phase. The 
growth of the purely growing drift instability would be an essential part of the 
proposed growth phase of the magnetospheric substorm. 

Some observational evidence supporting the auroral breakup mechanism 
proposed here is available. First, as noted above, a temperature gradient large 
enough to cause the growth of the purely growing electrostatic drift instability 
is likely to exist at times at the inner edge of the electron plasma sheet. The 
inner edge of the plasma sheet does respond to magnetic activity in a manner 
similar to the behavior of the auroral oval and the position of the inner edge is 
certainly closely related to the auroral zone. The location of the breakup of a 
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small isolated substorm has been placed at the inner edge of the plasma sheet by 
Johnstone, et. al. (1973). The existence of an electric field along the m^etic 
field lines has been cited by some investigators as consistent with observations 
of precipitated electron fluxes measured in breakup aurora (Johnstone, et. al., 
1973; Whalen and McDiarmid, 1972; Pongratz, 1972). However, the precipitated 
fluxes expected from the existence of a parallel electric field requires a detailed 
analysis and greatly depends on the shape and location of the electric potential 
along the field line and the manner in which the field is established. The precipi- 
tation of particles by suddenly turning on a parallel electric field has been investi- 
gated analytically by Chamberlain (1969). Although the mathematical device of 
suddenly turnii^ on a parallel electric field does not exactly correspond to the 
mechanism discussed here it is somewhat similar, A numerical investigation 
of the precipitated fluxes resulting from suddenly turning on a parallel electric 
field (after Chamberlain) has been performed by Matthews and Pongratz (1973). 
They use an initial energy spectrum similar to the plasma sheet spectrxun and 
a parallel electric field confined to the region near the magnetic equator. They 
are able to reproduce the precipitated fluxes measured using a sounding rocket 
in a breakup aurora. It might be expected that the mechanism proposed here 
should lead to similar precipitated spectra, but a detailed numerical calculation 
is necessary to support this assertion. 

The contention that the current limit is exceeded in the intense breakup arc 
may also have some observational support. Observations of precipitated particle 
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fluxes in a moderate substorm after the onset of breakup have been reported by 
Whalen and McDiarmid (1972). They report that the field-aligned current in the 
northern-most arc (the breakup arc after the onset of poleward expansion) ex- 
ceeds the critical current density given by Swift (1965) which is higher than the 
critical limit given by Kindel and Kennel (1971). They further observe that there 
is no correlation between proton and electron precipitation in the northern-most 
arc, a result v/hich might be expected from the mechanism proposed here. It is 
also interesting to note that Johnstone, et. al. (1973) claim the precipitated elec- 
tron fluxes observed spatially outside of the breakup arc are apparently due to 
"relatively strong" pitch angle diffusion. As pointed out above the growth of the 
drift instability should lead to enhanced pitch angle diffvision. 

It is appropriate to review some of the approximations used to analyze the 
growth of the electrostatic drift instability. The electric potential associated 
with the instability has been assumed to vary azimuthally as exp \ i and to be 
slowly varying in ct , The inner edge of the plasma sheet is known not to be 
azimuthally symmetric (Frank, 1971), The fastest growing mode, however, 
does have a large value and thus a short azimuthal wavelength. If the scale 
size of the azimuthal variations is much larger than the wavelength of the fastest 
growing mode then the assumption of azimuthal symmetry may be appropriate. 
Since the temperature gradient at the inner edge of the plasma sheet occurs 
over a narrow region (1-3 ) variations in the potential with a will probably 

be important. The solution for the electric potential in the curved field line 
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geometry is a formidable problem and would only merit consideration if the 
breakup mechanism proposed here appears to be correct. The approximation 
of electrostatic perturbation fields is not accurate near the magnetic equatorial 
plane. Before the growth of the drift instability a diamagnetic current exists at 
the inner edge of the plasma sheet due to the temperature gradient. The growth 
of the instability diminishes the diamagnetic cxirrent as exp-yt, which implies 
the existence of a magnetic field B’ with an associated electric field 

E' fvxE' = -1 

\ c 3t/ 

which the electrostatic approximation ignores. The magnitude of B’ is related 
to the diamagnetic current, 

and thus B’/B = 47Tp^ /B^ = where p^ is the electron pressure so that 
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but cE/B so E' ~ /3*E. Since /3* , the ratio of electron pressure to magnetic 
pressure is on the order of unity only near the magnetic equator, the magnetic 
field associated with the relaxation of the diamagnetic current will be important 
only near the equator. The growth of the electrostatic drift instability also has 
a large effect on the motion of particles mirroring near the equator. A thorov^h 
analysis of the purely growing drift should take into account the effects of the 
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magnetic field associated with the relaxation of the diamagnetic current on par- 
ticle motions. Hagege, et. al. (1973) have foimd in studying the "ion trapped 
mode" at higher frequencies eco) that finite / 3 * effects are stabilizing. 

However, /3* appears in integrals over the field line which tends to mitigate the 
problem since / 3 * decreases rapidly away from the magnetic equator. 

An important question concerns the ordering used to analyze the Vlasov 
equation. The ordering used here assumes a scale length of the same order of 
magnitude for the perpendicular electric field, the parallel electric field, and 
the temperature gradient. Frank (1971) asserts the pressure scale height is 
.4 R^. The fastest growing mode has a short azimuthal wavelength (-v.1 R ) anH 

c 

the parallel electric field will have as small a scale height as the ordering 
allows. Thus the ordering of scale heights used appears to be appropriate. The 
finite Lamor radius ordering should be used if the parallel electric field has a 
scale height much longer than the perpendicular scale heights. An investigation 
of the stability of the temperature gradient at the inner edge of the plasma sheet 
using the finite Lamor radius ordering has been carried out by Coroniti and 
Kennel (1970, see Chance et. al, (1973) for finite ft* corrections) in a straight 
field line geometry. They find an electromagnetic drift instability with a fre- 
quency on the order of 10“ i sec"i and a growth rate on the order of 10 “^ sec“^ 
having perpendicular wavelengths on the order of the ion gyroradius. This faster 
growing electromagnetic instability would, of course, relax the temperature 
gradient at the inner edge before the purely growing electrostatic drift instability 
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could grow. The electromagnetic drift instability may not occur in the midnight- 
evening sector of the magnetosphere, however. In fact, Coroniti and Kennel sug- 
gest the instability exists only in the post midnight sector (0200-1000 LT). The 
electromagnetic instability can be expected to grow when the temperature scale 
height is .4 only if the electron temperature is about 4 times the ion tempera- 
ture. Satellite measurements (Frank, 1967), at least for the evening sector, 
indicate proton temperatures of about 10-20 keV for the ring current and elec- 
tron temperatures near 1-5 keV (Vasyliunas, 1968; Schield and Frank, 1970). If 
the ion to electron temperature ratio is on the order of four to one, the scale 
height of the temperature gradient would have to be about 10"^ R^ for the electro- 
magnetic instability to grow. If the ion temperature is larger than the electron 
temperature and the scale height of the temperature gradient is near .5 R^ the 
electromagnetic drift instability investigated by Coroniti and Kennel cannot be 
expected to grow whereas the purely growing electrostatic drift instability can 
grow. 

The growth rate of the purely growing drift instability has been estimated 
as 7 ~ -tiDp, where ^ is small enough that . Since magnetospheric ions 

have a higher temperature than electrons I can be on the order of 10 ^ without 
violating the second adiabatic invariant, J, of the ions and y ^ 10"^ sec"^. The 
loss of electrons to the atmosphere modifies the growth rate but this effect is 
ignored here. Kennel (1969) estimates the electron life-time is less than 10^ 
seconds at the inner edge of the plasma sheet during strong diffusion. Strong 
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diffusion must be assumed absent if the purely growing drift instability is to grow 
although as has been previously noted the growth of the instability should sub- 
sequently lead to strong diffusion. The drift instability can have a growth rate 
on the order of (10 minutes)"* . This growth rate may be fast enough for the 
purely growing drift instability to near saturation during the growth phase of 
the magnetospheric substorm. McPherron (1972) gives the duration of the growth 
phase as about 1 hour. If the parallel electric field of the instability at saturation, 
E]| g = <l)g/L, is estimated by 

''ll 

the potential amplitude along the field line near the equator before breakup may 
be on the order of kT/e, a few kilovolts. 

In summary, the mechanism for auroral breakup proposed here depends on 
a sequence of related events. The temperature gradient which exists at the inner 
edge of the electron plasma sheet is the source of energy for a purely growing 
drift instability (the pressure gradient appears to be strong enough to meet the 
requirements for instability). Initially, precipitation of electrons by whistler 
noise in the vicinity of the inner edge is assumed weak. The drift instability 
grows (sometime during the growth phase of the magnetospheric substorm) with 
little particle precipitation, alters the equatorial velocity space distribution 
creating a parallel electric field which is concentrated near the magnetic equator. 
The alteration of the equatorial velocity space distribution is an anisotropy in 
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pitch, angle wiiich leads to generation of whistler noise and pitch angle diffusion 
of electrons. The growth rates and band width of the whistler noise increases 
as the parallel electric field of the drift instability gets stronger. The drift in- 
stability continues to grow until the electron loss by whistler generated pitch 
angle diffusion is large. If the current in the ionosphere due to electron loss is 
strong enough, the current limit is exceeded, and turbulent resistivity set up in 
the ionosphere leading to rapid and intense electron precipitation. This rapid 
electron precipitation when the current limit is exceeded is proposed to be the 
auroral breakup. If the current limit is not exceeded the parallel electric field 
is relaxed by strong pitch angle diffusion, and no breakup occurs. 
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Figure Captions 


Figure 1. Directional differential energy spectrums of electron intensities 

measured by OGO 3 in the midnight sector of the magnetosphere near 
the equatorial plane (Schield and Frank (1970) figures 4 and 7). The 
figures show the changes in the spectrum observed on satellite passes 
across the inner edge of the plasma sheet. 

Figure 2. Energy densities of both protons and electrons for the energy range 

1-50 keV as measured by OGO 3 (Frank (1971) figures 1 and 5), show- 
ing the location of the trapping boimdary and plasmapause. Figure 2a 
corresponds to the detailed data shown in figure 3. 

Figure 3, Directional differential electron intensities for various energy band- 
passes as functions of L value. The center energies of the bandpasses 
are noted in the figure. Note the increase in low energy particles for 
decreasing L value as the inner edge of the plasma sheet is crossed. 
The data for this figure is taken from Shield and Frsink (1970) figure 2. 

Figure 4. Directional differential energy spectrums of electron intensities cor- 
responding to figure 3 at 6.6 and 6.8 . The two curves are the 

intensities for a Maxwellian energy distribution at the temperatures 
shown (1.0 and 1.55 keV) for a constant number density. Note the 
lack of low energy particles for the measured plasma sheet intensities 
compared to the intensify for the corresponding Maxwellian. 
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